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We consider the low energy collective mode spectrum of 
a superfluid Fermi gas in a spherical trap in the collisionless 
regime. Using a self-consistent random-phase approximation, 
the effects of superfluidity on modes of dipole and quadrupole 
symmetries are systematically examined. The spectrum is 
calculated for varying pairing strength and temperature and 
we identify several spectral features such as the emergence 
of Goldstone modes that can be used to detect the onset of 
superfluidity. Our analysis is relevant for present experiments 
aimed at observing a superfluid phase transition in trapped 
Fermi gases. 

Pacs Numbers: 32.80.Pj, 05.30.Fk, 67.55. Jd 

Properties of trapped Fermi gases are attracting in- 
creasing attention in the field of ultracold atomic sys- 
tems. Experimentally, the trapping and cooling of 
fermionic alkalis has been demonstrated reaching tem- 
peratures as low as - Ti./4 for ^Or § and ^Li |-|] 
with Tp denoting the Fermi temperature. One goal of 
the present experimental effort is to observe a transition 
to a superfluid phase at some critical temperature Tc for 
systems with an attractive interaction between atoms in 
two different hyperfine states [g|. Several proposals for 
the detection of this phase transition using the interac- 
tion with light Ig , thermodynamic and various collective 
and quasi particle properties [Q have been presented. 

The purpose of this paper is to examine the effects 
of superfluidity on the low energy collective mode spec- 
trum. The modes can be measured with high precision 
as shape oscillations of the atomic cloud and such studies 
have proven to be very useful in the field of Bose-Einstein 
condensation ^. Ultracold atomic Fermi gases provide 
an experimental realization of a neutral superfluid with 
a simple tunable interaction. Both regimes where the 
pairing is a small perturbation on the normal properties 
of the gas and where it dominates the low frequency re- 
sponse could be experimentally obtainable. Also, finite 
temperature effects should be measurable. This makes 
the study of collective modes in these systems interesting 
as they in some sense are simpler and easier to experimen- 
tally manipulate than other well known superfluid Fermi 
systems such as liquid ''He and the atomic nucleus. 

We present an analysis of the low energy collective 
mode spectrum of a trapped superfluid Fermi gas for var- 
ious temperatures and coupling strengths. The trapping 
potential has a qualitative effect on the low energy modes 
as opposed to the high energy regime, where one can use 
the well-known results for a homogeneous superconduc- 
tor 0. We assume that the underlying normal gas is 
in the collisionless limit such that the lifetime r of the 
quasi particles is much longer than the characteristic pe- 
riod of motion (i.e. lott 3> 1 for atoms in a harmonic 



trap of frequency lot)- It is shown that superfluidity has 
qualitative effects on the modes of dipole and quadrupole 
symmetries and we thus propose several ways of detect- 
ing the superfluid transition. An essential difference be- 
tween the density and spin-density fluctuation modes is 
demonstrated and we also analyze the effects of a non- 
zero temperature. 

Consider a gas of fermionic atoms of mass m in two 
hyperfine states \a =t, i) trapped by a spherically sym- 
metric harmonic potential t^o(r) = imuj^r'^. The num- 
bers Na of atoms trapped in each hyperfine state are 
assumed to be equal as this is the optimum situation for 
Cooper pairing. The quasi particle excitations i] of the 
system are calculated by the Hartree-Fock-Bogoliubov 
(HFB) method. Using a zero-range pseudo-potential 
method to describe the interaction between the atoms, 
the Bogoliubov-de Gennes (BdG) equations become pOJ : 



E„u^{r) = [Ho + Wir)u„ir) + A(r)w^(r) 
E^,v^{r) = -[no + W{r)]v^ir) + A(r)w^(r). 



(1) 



Here, Hq = ~{h^/2m)V^ + Uo{r) - fip, W{r) = 
.g(V'J(r)'0o-(r)) = diPai^)) is the Hartree potential, and 
ijja-i^) is the field operator that annihilates an atom 
in spin state a at position r. The coupling constant 
is g = Airah^ /m with a being the s-wave scattering 
length between atoms in the two different hyperfine states 
(a < 0, i.e. attraction) and ^ip denotes the chemi- 
cal potential. The pairing field is defined as A(R) = 
-5lim^^o5r.[r(Vi|(R-Hr/2)i/);(R-r/2))] fill. The quasi 
particles with energies i?,, are described by the Bogoli- 
ubov wave functions M,,(r) and w,;(r). 

We are in this paper interested in the response of the 
gas in the superfiuid phase to various particle conserv- 
ing time-dependent external perturbations of the kind 
F[t) ex exp(ia;i) Yla I d^'rFcr{r)pcr{v). In the collisionless 
limit, the appropriate framework to calculate the collec- 
tive mode spectrum is the self-consistent random phase 
approximation generalized to superfluid systems |12| , |l3| . 
The linear response of the superfluid is characterized by 
a matrix consisting of two-particle correlation functions: 



H(^) 



((PtPt)) ((PTPi» 

((PiPT» ((PiPi)) 

((XPT» ((xPi» 

((X^PT» ((Xt/5i» 



((PTX)) ((PTX^)) 

{{PIX)) {{Pim 

iixx)) iixx^)) 

iix^x)) iix^x^)) 



(2) 



Here, {{AB)) is the Fourier transform of the retarded 

fimction ~iQ{t - t'){[A(Y, t), B{v' , t')]) and (...) denotes 

the thermal average. The operator x ~ ^i(r)^T('^) de- 
scribes the local fluctuations in the pairing fleld. Note 



that we have spht the correlation functions into the differ- 
ent hyperfine components since we also wish to consider 
spin modes in which the two hyperfine species oscillate 
in anti-phase. For the normal phase, the density and 
spin-density fluctuation modes can be obtained by calcu- 
lating the upper left quadrant ofIl{Lu) in Eq.(||) [Q. For 
the superfluid phase however, fluctuations in the pairing 
field are coupled to density fiuctuation modes through 
the non-particle-conserving correlation functions such as 
{{px)) Ei-nd we need to consider the full matrix in Eq.(y). 
Within RPA, we have 



u{oo) = [i-no{uj)g]-^Uo{cu) 



with 



|j(r-r') 



10 
10 
1 
10 



(3) 



(4) 



describing the interaction between quasi particles in op- 
posite hyperfine states. Here the matrix products denote: 
AB =J(fr"A{r,r")B{r",r') and Uo{uj) is the matrix 
in Eq. (g) calculated for non-interacting excitations using 
the HFB approximation given by Eq.(n|). 

For a spherically symmetric trap, the correlation func- 
tions split into terms describing the various multipole 
modes ||l^]. For instance, we obtain for the correlation 
fimction ((p^(r)/5^(r'))): 



((P.(r)p.(r')» 



LM 



EL{r,r\u;)YLM{0,4>)Y£j,j{e',cj)') 



(5) 



with r denoting the distance from the center of the trap 
and the Ylm{0, </>) the usual spherical harmonics. The ex- 
pressions for the other correlation functions appearing in 
Eq.(ph are completely equivalent. Using Wicks theorem, 
we express the independent particle correlation functions 
((. . .))o used to form no(a;) in terms the quasi particle 
energies/wave functions obtained from Eq.(n|). For in- 
stance, for {{p„{r)pc{r')))o, we obtain 



^OL{'r,r' ,uj) 



IV 



E 

rrn' 



47r(2L+ l)r'r 

uuuu{f — /') uuvv{f + /' — 1) 

^ E + E' -hw-iS 

vvvv{f' - /) 



E-E' -huj-iS 

vvuu{l - f - f) 

' E + E' + hLU + i5 - 



-E + E' -huj- iS 



(6) 



Here, we write the Bogoliubov wave functions from 
Eq.(|l|) with angular momentum quantum numbers Im 
and energy £",,/ in the form «,,;,„ (r) = r"iu,,;(r)Yim(6', (j)) 
and Vriimi'^) — r^^v,-iiir)Yim{0,<P)- We define uuvv = 
Uni{r)uni{r')vr,'i' {r)vn'i' {r') etc., and / = [exp{/3Er,i) + 
1]~^ is the Fermi function with (3 — l/ksT. The Clebsch- 
Gordan coefficients {ll'00\LO) yield strong selection rules 



for / and I': for the monopole modes {L — 0) one has 
I ~ /', for the dipole mode (L = 1) /' = Z ± 1 etc. 

There is one last technical issue to resolve: The correla- 
tion function ((x(r)x^(r'))) is ultraviolet divergent when 
expressed as a sum over Bogoliubov wave functions as in 
Eq.(6) [1^]. Using the pseudo-potential method in a way 
similar to that described in Ref. |ll], this divergence is 
removed by the substitution 



((X(r)xt(r'))> -> {{x{r)xHr'))) + x<5(r - v')GZ:{r). 



(7) 



Here G!," 



r) is the part of the the single particle Greens 
function G^j^(r,x) = {r + x/2\HQ^\r + x/2) that diverges 
as l/x for a; ^ |1^. One can show that for the homo- 
geneous case, this procedure is equivalent to subtracting 
the non-regularized gap equation from the calculation of 
iixX^)) Etnd the diagonalization of n(ci;) then yields the 
well-known Bogoliubov- Anderson mode ||l^ . 

The structure of the self-consistent RPA calculation 
is then the usual one: First we obtain a self-consistent 
solution to the BdG equations given by Eq.(|l|); then 
Ho(a;) is formed using the obtained Bogoliubov wave 
functions/energies; and finally the RPA response func- 
tion is calculated from Eq. (0). The poles of H(u;) then 
give the collective modes of the gas. 

The quantity of experimental interest is the strength 
function 



S{F,co) = J2 



^-0E„ _ g-/3£„ 



\{n\F\m)\^S{huj 



En — Em) 

(8) 



which is directly related to the net transitions per unit 
time with energy fko induced by the operator F. Here Z 
is the grand partition function and \n) is an eigenstate 
of the Hamiltonian with energy En . For operators of the 
form F{t) oc exp{iujt) ^^ / d^rF^{r)pa{r), it is given by 



S{F,u:) = -^Y.l d\<fr'F,F'^,Im[{{p,p',,))]. 



(9) 



We shall calculate it for ^^(r) = F^(r)YLM{0,i^) with 
F|(r) ~ Fi{r) = r^ for the quadrupole {L = 2) den- 
sity mode. For the dipole symmetry (L — 1), we take 
-^t(''") ~ ~-^i(''') = '^ exciting the lowest spin-dipole 
mode. Taking F|(r) = Fj(r) = r for i = 1 would sim- 
ply excite the center-of-mass mode at w = ut- Taking 
L = excites modes of monopole symmetry. The infiu- 
ence of superfluidity on the monopole modes is somewhat 
complicated and will be analyzed in a future publication. 
The properties of the superfluid gas in general depends 
on quantities such as the size of the coherence length f of 
the Cooper pairs, the quasi particle spectrum obtained 
by solving Eq. (p, Hlot and the mean field splitting of 
the trap levels. For a large number of particles trapped 
and a strong interaction one has A(r) 3> ficvT for small 
r and therefore ^(r) ^ R with R denoting the spatial 



extent of the cloud. As we shall demonstrate, the pair- 
ing has strong effects on the low energy spectrum in this 
regime and should thus be relatively easy to detect. We 
therefore believe this limit to be of most relevance for 
the present experimental effort to observe the superfluid 
transition. To illustrate some important physical con- 
cepts, we also present results for weaker coupling where 
the pairing mainly occurs in the levels of the oscillator 
shell that is closest to the chemical potential (valence 
shell). In this case, there are quasi particle excitations 
with Ejj^Qi < Hlot since the pairing energy is smaller 
than hujT- However, we do not in this paper consider the 
theoretically interesting A <^ Hcot limit in detail with A 
denoting the pairing energy. In this regime the interplay 
between the shell structure of the trap and the Cooper 
pairing leads to a more complex phase diagram. 

We first present typical results for a system with rel- 
atively weak pairing. Figure (0) exhibits S{F, oj) (in 
trap units) given by Eq. (|9|) for the spin-dipole, and 
quadrupole modes. There are ~ 1.4 x 10'* atoms trapped 
with a coupling strength such that the pairing field in 
the center of the trap is A(r = 0) ~ 0.69hu}T for T = 0. 
For the specific parameters chosen, it is mainly the levels 
within the harmonic oscillator shell with principal quan- 
tum number n = 33 [with energy £,n = {n + 3/2)hLUT in 
the non-interacting case] which Cooper pair. The Hartree 
field has lowered the energy of this shell so that it is sit- 
uated around the chemical potential at fj,p = 32hujT- 
We have added a small imaginary part F = O.OOSfiwT to 
the frequency to model a smooth response typically ob- 
served experimentally. For comparison we also plot the 
response for fc^T = Hcot where the gas is in the normal 
phase {kBTc ~ 0.9/iwt) and both the quasi particle and 
the collective (RPA) response is shown. As expected, 
there is a large response at a; ~ ujt,2ojt for the spin- 
dipole and quadrupole response respectively for the gas 
in the normal phase. The shift of these peaks away from 
the ideal gas result is, of course, a consequence of the 
particle-particle interaction |14| ; we will not discus these 
shifts further in this paper as we are concentrating on 
features specific to superfluidity. 

From Fig. (|l|) , we see that the presence of superfluidity 
introduces a qualitative new feature in the quadrupole 
response: a low frequency response for lo of the order 
'^ OAujt- The low energy response essentially comes from 
the breaking of Cooper pairs in the lowest quasi particle 
shell. As pointed out by Baranov |1^], the low energy 
quasi particle excitations [£",,=0,; ^ A(r = 0)] are surface 
modes since the pair correlations imply that the interior 
of the paired system is an incompressible medium which 
expels the low energy quasi particles. Reflecting this ef- 
fect, the response observed at ~ OAHujt corresponds to 
a pair of quasi particles moving in the outer part of the 
cloud; the strength of this response in RPA is increased 
and the frequency decreased as a result of the (weak) 
coupling with the collective mode at '^ 2hujT- This is 
illustrated by the inset in Fig. (^ which compares in 
detail the low energy single particle response calculated 



from no(w) with the collective response. The single par- 
ticle response is generated by breaking of Cooper pairs 
due to excitations of the kind jl^Qit'jL^Qii A^o) with 
{I' = l,l±2) costing E,,=oi + E^=oi, in energy. Here, |$o) 

is the superfluid ground state and -yLg. creates a quasi 
particle with quantum numbers rjla. The excitation en- 
ergy associated with the lowest quasi particle states can 
mainly be attributed to the pair breaking energy A and 
we write Eri=oi + Ej^^qii ^ 2A. Thus, the low energy re- 
sponse consists of pair breaking excitations with lu '^ 2A 
which is ~ OAHlut for this set of parameters. As we 
shall see, this response is the precursor of a Goldstone 
mode present at stronger pairing. Furthermore, the peak 
at ~ 2ujt is fragmented and partially damped due to 
the pairing. We will demonstrate below that this mode 
simply disappears with the emergence of the Goldstone 
modes for stronger pairing. 

For the spin-dipole mode, we see that the main effect 
of superfluidity in the A < Hlot regime is that it frag- 
ments and partially damps the normal phase response at 
~ ujt- As for the quadrupole mode at ^ 2ujt, this ef- 
fect is the precurser of the situation in the strong pairing 
regime, where we will show that there is no well-defined 
low frequency spin-dipole mode. 

Detailed effects of superfluidity on the mode spectrum 
in the weak pairing regime, such as the way the normal 
phase quadrupole/spin-dipole modes at ~ 2wt/'^t are 
fragmented and damped with increasing pairing, depend 
sensitively on the specific parameters of the system such 
as the exact location of the chemical potential relative to 
the quasi particle bands. We will discuss these effects in 
detail in a future publication. 

We now consider the stronger coupling limit when 
^ < i?. Figure (|) exhibits S{F, uj) for - 1.6 x 10* atoms 
trapped with a coupling strength such that the pairing 
field in the center of the trap is A(r = 0) = 6hujT for 
T — 0. We have taken F — O.OObhujT and the response 
is calculated for T = and ksT = IAHujt where the gas 
is still superfluid and ksT = S.OHujt ( where the gas is 
in the normal phase [UbTc — 2.8hujT)- Again, we plot 
both the quasi particle and the collective response. A 
semi classical analysis shows that the lowest energy quasi 
particles are localized near the minimum of the function 
A{r)^ + [Uo{r) + Wir) + tf{l + l/2)y{2mr^)-HF]^ 0. 
With increasing pairing, the width of this minimum nar- 
rows and the lowest quasi particle energy can be shown to 
increase as ln{fip/huJT)- For the specific parameters cho- 
sen above, the lowest quasi particle states have energies 
En=oi ^ hujT- 

From Fig. (||), we see that for the quadrupole mode 
there is for T — now only one low frequency mode sit- 
uated at w ~ IAujt- This lower mode comes from the 
pairing degree of freedom [the lower right 2x2 part of 



the matrix H in Eq.( 
known Goldstone moc 



J)]; it is an example of the well- 
es [|T^ with a finite frequency due 
to the trapping potential. The low energy quasi parti- 
cle response observed for weaker pairing [Fig. M)] has 



now become a collective (Goldstone) mode. For these 
modes, one can develop a hydrodynamic theory valid for 
^ <C -R IQ predicting a quadrupole Goldstone mode at 
ix> = \/2ujT fl^l in agreement with our result. Of exper- 
imental relevance is the temperature dependence of the 
spectrum. From Fig. (||), we see that as T is lowered 
through Tc the mode at lj ~ 2u!t is damped and a low 
frequency response emerges. At intermediate tempera- 
tures <T <Tc, the collective mode may suffer appre- 
ciably Landau damping [see the spectrum for T = Tc/2 
in Fig. (0)]. The spectrum in Fig. ph is consistent with 
an interpretation in which the collective mode emerges 
as an undamped excitation as soon as the temperature 
is low enough so that the lowest quasi particle excitation 
energy is > V^Hiot- For the parameters in Fig. (0), this 
occurs for T - Tc/A. 

We finally consider the low lu spin-dipole response. 
From Fig. (0), we see that the normal phase mode at 
~ LOT is totally damped in the superfluid phase. This is 
to be expected for A > hajr- The only kind of motion 
possible for the superfluid is potential flow, i.e. density 
fluctuation modes jlq]. All other modes such as spin- 
density modes requires the presence of a normal compo- 
nent which decreases with T — > 0. The superfluid cannot 
participate in a mode of relative motion of the two hyper- 
flne states; one needs to break pairs and the response lies 
above the pairing energy. The damping of the spin-dipole 
mode can also be understood by writing the spin-dipole 
operator in terms of the quasi particle operators: [p-\ (r) — 

P4(r)]|$o) oc E„r,'K(r)wv(r) - Wr,(r)ur,'(r)]7^^7j,J$o). 
Pairing tends to mix a particle state with its time re- 
versed hole-state yielding Wr)(r) ~ w^(r) for a low energy 
excitation and therefore [u,,(r)w^/(r) — w,,(r)u,,/(r)] ~ 0. 
Thus, the low lo spin-dipole mode excited by an odd op- 
erator under time-reversal will gradually disappear for 
T ^ as can be seen from Fig. ^. 

In conclusion, we have identified and discussed the 
physical origin of several clear-cut spectral features of su- 
perfluidity on the low frequency collective mode spectrum 
of dipole and quadrupole symmetry in the coUisionless 
regime for weak and strong pairing. For a weak pairing 
A < hujT , the main effect of superfluidity is the presence 
of a low Lo response for the quadrupole symmetry. Also, 
the normal phase quadrupole and spin-dipole modes at 
~ 2ujt and ~ lut respectively are partially damped. For 
stronger pairing with ^ <^ R, the Goldstone modes dom- 
inate the low frequency spectrum and the quadrupole 
mode is shifted to w \/2wt in agreement with hydrody- 
namic theory. In this regime there is no well-defined low 
energy spin-dipole mode. As illustrated above, the emer- 
gence of these effects should be straightforward to detect 
as one lowers T through T^ and we believe the results pre- 
sented to be relevant for the present experimental effort 
related to observing the predicted BCS transition. 

We would like to acknowledge valuable discussions 
with C. J. Pethick. 
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FIG. 1. The collective (RPA) and quasi particle (Q.P.) re- 
sponse for the spin-dipole and quadrupole modes for A < Itut 
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FIG. 2. The collective (RPA) and quasi particle (Q.P.) 
response for the spin-dipole and quadrupole modes for 
A > fkJT- The insets show various frequency regions in detail. 



